Running couplings for the simultaneous 
decoupling of heavy quarks 



Steven D. Bass*' 1 , R.J. Crewther b , F.M. Steffens c ' d ' 2 , 

A.W. Thomas 6 ' 3 

a Institute for Theoretical Physics, Universitdt Innsbruck, Technikerstrasse 25, 

A 6020 Innsbruck, Austria 

h 'Department of Physics and Centre for the Subatomic Structure of Matter 
( CSSM), University of Adelaide, SA 5005, Australia 

c Instituto de Fisica Teorica, UNESP, Rua Pamplona 145, 01405-900 Sao Paulo, 

SP, Brazil 

d Mackenzie Presbiteriana University, FCBEE, Rua da Consolacao 930, 01302-907 

Sao Paulo, SP, Brazil 

, Jefferson National Laboratory, 12000 Jefferson Avenue, Newport News, 

VA 23606, U.S.A. 



Abstract 

Scale-invariant running couplings are constructed for several quarks being decoupled 
together, without reference to intermediate thresholds. Large-momentum scales can 
also be included. The result is a multi-scale generalization of the renormalization 
group applicable to any order. Inconsistencies in the usual decoupling procedure 
with a single running coupling can then be avoided, e.g. when cancelling anomalous 
corrections from t, b quarks to the axial charge of the proton. 

Key words: Renormalization group, thresholds, mass logarithms, matching 

conditions, axial anomaly 

PACS: ll.10.Hi, 12.38.Aw, 12.38.Cy, 12.39.Hg 



Normally, heavy quarks are decoupled one at a time, in order of decreasing 

mass m q , q = t,b, First the top quark is removed, with the b- and c-quark 

thresholds left intact. Then b is decoupled without destroying the c-quark 

1 Supported by the Austrian Science Fund, contract P17778-N08. 

2 Supported by FAPESP (03/10754-0) and CNPq (308932/2003-0). 

3 Supported by DOE contract DE-AC05-84ER40150, under which SURA operates 
Jefferson National Laboratory. 



Preprint submitted to Physics Letters B 



2 February 2008 



threshold, and so on: 



m t — > oo first, then nib — > oo, then m c — > oo. (1) 
This sequential approach has two basic problems: 

(1) On a logarithmic scale, m t may be sufficiently large compared with m^, 
but not 4 rrib compared with m c . 

(2) More seriously, it is inconsistent to retain both 0(d(m t )) and 0(a(mb, c )) 
corrections when the first step m t — > oo assumes a(m t ) <C a(mj iC ); here 
a(M) is the usual running coupling at scale M. 

Our proposal is to replace the standard procedure by a new technique based 
on an extended renormalization group (RG) where heavy-quark decoupling 
is genuinely simultaneous and all mass logarithms are summed together. We 
show the need for this with an explicit example: anomaly cancellation for t, b 
quarks decoupled from the axial charge of the proton. 

We let two or more mass logarithms grow large together relative to a light- 
quark scale p,, e.g. 

\n{m c /p) < \n(m b /ji) < \n(m t /fi) -> large (2) 

for three heavy quarks. These heavy quarks are decoupled together, i.e. in 
a single step, without reference to intermediate thresholds. Instead of a single 
running coupling a(M), there are separate running couplings a t , a b , ■ ■ ., one for 
each heavy quark h = t,b, .... The leading asymptotic power of an amplitude 
A is expanded as a multinomial in a t , aj,, 

A = 52at<4...Au... (3) 
kt.. 

with all {rrih} dependence contained in a t ,ai,,.... The coefficients Am... are 
amplitudes produced by the residual light-quark theory. They can be identified 
by observing that the special case 

\a(m t /jl) > \a(m b /jx) > . . . (4) 

of the simultaneous limit (2) must give the same results as a direct application 
of the sequential limit (1). This technique depends on the fact that the residual 
theory is the same, irrespective of how the heavy quarks are decoupled. 

4 The t quark, the W ± and Z° bosons, and (presumably) the Higgs boson, provide 
another example of this problem. When more than one of these particles sets the 
scale, or there are momenta of similar magnitude, it is not reasonable to suppose 
that one of the corresponding logarithms dominates the others. See [1,2,3,4,5] for 
alternative methods. 
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A key step is the construction of simultaneous running couplings at, a b , . . . 
for the general limit (2). Each ah is RG invariant, depends on all of the 
mass logarithms in Eq. (2), and in the limit (4), can be expanded in running 
couplings a t , a b , and (if needed) a c associated with separate steps in the 
sequential limit (1). 

Let af = <?j/47r and m q f be the strong coupling and quark masses renormalized 
in the presence of / flavours and three colours according to the MS (modified 
minimal subtraction) scheme. We need the functions 

(3 f (x) = ~(33 - 2f) - ^(153 - 19/) + 0(V) , (5) 

S f (x) = + 0(x 2 ) (6) 

associated with the Callan-Symanzik operator 

where p, = VAne^^^im is the MS scale (//dim being the scale for dimensional 
renormalization). Note that p, is chosen to be independent of /, that it sets 
the fixed low-energy scale needed to define simultaneous limits, as in Eq. (2), 
and that the heavy-quark masses for these limits are renormalized in the six- 
flavour theory: 

m h = m m , h = t,b,... (8) 

The running coupling a t for the heaviest quark t is defined by adapting Wit- 
ten's prescription [6] for a single heavy quark. We set / = 6, 



rat 

\n{m t /p) = / dx{l-5 6 (x)}/(3 6 (x) (9) 



and require that a 6 be matched directly to the coupling af of the residual 
f-flavour theory in the limit (2). By construction, a t is scale invariant with 
respect to the original / = 6 theory: 

V 6 a t = . (10) 



Scale-invariant definitions of the other running couplings are then easily pro- 
duced. The idea is to exploit a characteristic feature of simultaneous decou- 
pling — the presence of large scale-invariant logarithms such as 5 

\n(m t /m b ) = \n(m t /p) - \n(m b /p) . (11) 

5 In the standard procedure (1), one chooses a scale for each step such that large 
logarithms never appear [7]. For simultaneous limits (2), that is no longer possible. 



3 



This logarithm can be used to set the relative scale between a t and a b , 

ln(m t /m b ) 6 = £*dx{\ - S 5 (x)}/p 5 (x) , (12) 

with a c defined similarly if the c quark is also being decoupled: 6 

In (m b /m^j g = J dx {1 — 8i{x)} / f3±(x) . (13) 

Note that a b and a c exhibit the full / = 6 scale invariance of the original 
theory 

V 6 a b = = V 6 a c (14) 
irrespective of the integrands chosen in Eqs. (12) and (13). 

We will show that these definitions produce simultaneous running couplings 
a t , a b , . . . with the property 

a h = a h + 0(al K ) (15) 

in the sequential limit (4), where a t , a&, • • ■ are Witten's running couplings for 
successive steps in Eq. (1): 

\n(m t6 /p) = [ dx{\ - 5 6 (x)}/(3 6 (x) , (16) 

J Old 

\n(m b5 /fl) = [""dx {1 - S 5 (x)}/p 5 (x) , (17) 
ln(m c4 /At) = £ "dx {1 - 5 4 (x)}/fa(x) . (18) 

Eq. (15) says that at and ah coincide asymptotically in leading order (LO) 
and next-to-leading order (NLO) 7 . LO agreement ensures that any amplitude 
A which can be expanded in the sequential limit is analytic in {a^} when ex- 
panded in the simultaneous limit, as foreshadowed in Eq. (3), while agreement 
to NLO accuracy is convenient for practical purposes. 

For example, let us decouple t, b simultaneously from an amplitude A. We 
match the simultaneous expansion (3) of A to terms produced by sequential 
decoupling, where first the top quark is decoupled with m b held fixed, 

A = A + a t Ai + a 2 t A 2 + 0{a\) (19) 

6 These definitions make sense for all mass values consistent with m t ^ m c . 
This includes equal-mass cases, where running couplings coincide, e.g. at = a b for 
m t = m b . 

7 For at and at, the agreement is to all orders in the leading power: the defining 
equations (9) and (16) coincide. 
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and then b is decoupled from each of the residual five-flavour coefficient am- 
plitudes Ak, k — 0, 1, 2, . . . to produce four-flavour amplitudes Am'- 

A k = A k0 + a b Aki + a 2 b A k2 + 0(a 3 b ) . (20) 
Combine Eqs. (19) and (20), collecting terms of the same power in a t , b ' 
A = Aqo + {a t ^lio + Moi} + {a 2 A 20 + a t a b A u + a 2 b A 02 ) + 0(&l b ) . (21) 

Then, for any fixed value of ln(m t / jl) / ln(m b / jl) from 1 to infinity, the linear 
(LO) and quadratic (NLO) corrections for simultaneous decoupling involve 
the same coefficient amplitudes as in Eq. (21), because of the property (15): 

A s = Aqo + {a t A w + a b A 01 } + {a 2 t A 20 + a t a b A u + a 2 b A 02 } + 0(a 3 6 ) . (22) 

Coefficient amplitudes for higher orders can be calculated directly from cubic 
and higher-order corrections to Eqs. (15) and (21). Phenomenology based on 
single-limit expansions like (22) is perfectly consistent because all running 
couplings are of the same order, e.g. a t ~ a b - 

The proof of Eq. (15) depends on RG equations which relate the initial F- 
flavour and residual /-flavour theories. These equations are equivalent to the 
rule [8] T>p = T>j when these operators act on /-flavour quantities. Applied 
to «/ and mm q f, this rule yields Witten's relations [6], generalized to include 
simultaneous cases (2): 

Vpctf = fifiptf) , Vpinrnqf = Sf(otf) for all / < F. (23) 

An immediate consequence of Eq. (23) is that the sequential running couplings 
defined by Eqs. (16)— (18) are all RGi? =6 invariant [8]: 

V 6 a t = V 6 a b = V 6 a c = . (24) 



The proof also depends on matching conditions [8,9,10,11,12] for the / = 4, 5, 6 
couplings and masses in Eqs. (16)— (18). Let the decoupling of h correspond 
to reducing the number of flavours from / + 1 to /. Then a scale-invariant 
matching function Tf+\^f can be defined by the condition [8] 

f&h f<*h 

/ alx/(3 f+1 {x)= dx/P f (x)+F f+1 ^ f (a h ) . (25) 

Ja f+1 Ja f 

It contributes at NNLO (next-to-NLO): 

*> + w(fi*) = - 12^(33-2/) ^ + °(^) • (26) 



5 



There is also an invariant mass-matching function Gf+i^f, 

. m a (f + D r a h 5 fix) [ & h , 8f+i(x) _ , _ . 



which contributes at NNNLO (next-to-NNLO) [8]: 

89 
"432tt 2 



89 / \ 

0/+w(**) = -7ooZ2«fc + °(0 • ( 28 ) 



Now let Eqs. (12), (16), (17) and the identity 

ln(m t /TO 6 ) 6 = ln(m t6 ///) - ln(m 65 //i) - ln(m 66 /m 65 ) (29) 

be combined with Eqs. (25) and (27) for the case / = 5, h = t. The result is 
a relation between the various running couplings for t and b: 

[°t l-8 5 (x) p l-6 5 (x) 
The b, c running couplings are similarly related: 

/ dx a l \ = / dx a l \ V ^5-^4{a b ) ~ Gb-+A{Oi b ) . (31) 

Eq. (15) is an immediate consequence of Eqs. (30) and (31): the matching 
functions do not contribute at NLO, and we have7° at = a t for the top 

seq 

couplings. Indeed, the NNLO corrections are easily found: 

a b = a b - -^alat + 0(a t 4 b ) , (32) 



seq 72tt2 

^r^«c(«b + «*) + o(®t,b,c) ■ ( 33 ) 



a, = a. 



seq 72tt2 



As noted below Eq. (9), a 6 and m q6 have to be matched directly to their 
counterparts in the residual f-flavour theory. This is done in two steps: we 
replace Witten's invariant mass [6] rhh for the heavy quark h in Eq. (25) 

ln(m h /jl) = / dx/f3 f+1 (x) (34) 

Ja f+1 

with a set of simultaneous RGf=6 invariant masses m t ,m b , . . ., and then define 
simultaneous matching functions J-q^ and £/6->f- 

In the case of simultaneous t, b decoupling, where the residual theory has four 
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flavours, m t and mj are defined by the equations 

rat 

\n(m t /p) = / dx/[3 6 (x) , 

Ja G 
rat 

\n(m t /m b ) = / dx/(5 b {x) . (35) 

Ja b 

Then T§-, 4 is defined by the condition 

rab 

ki(m b /p,) = dx/p 4 (x) + Fe^at, a b ) . (36) 

J Q4 

Being RGp =6 invariant, ^6->4 can depend only on a t and a b . The matching 
condition between a 6 and «4 for simultaneous t, b decoupling is the result of 
eliminating a t and a b from Eqs. (35) and (36). 

An alternative definition of T%^ 4 is obtained by eliminating rn t and m b from 
Eqs. (35) and (36): 

rab 1 fat \ ra& \ 

Fo^at, a b ) = - / dx——- - / dx—— - / dx—— . (37) 

Comparison with Eq. (27) then suggests the following RGi? =6 invariant defi- 
nition of a simultaneous mass-matching function Qq^: 

In =/ dx——+ dx——+ dx—— + y 6 ^ 4 {a t , a b ) . (38) 

m£4 Ja 4 P±{X) Ja b (J 5 (X) J a t I3q(x) 

Eqs. (37) and (38) can be matched to Eqs. (25)-(28) and (32) in the sequential 
limit (4) to determine the leading terms in (NNLO) and Q§^ 4 (NNNLO): 

Fa-** = ~^7j^( ab + a *) + °( a tt) > 
89 / \ 11 

Q(i ^ 4 = ~432^V^ 2 + °$) ~ m^ abat + °( a *J ■ (39) 



For simultaneous t, b, c decoupling, we add 

ra b 

\n(m b /m c ) = / dx/(3 4 (x) (40) 

J a c 

to the definitions (35) and specify matching functions ^6-»3 (NNLO) and <?6^3 
(NNNLO) as follows: 

, m c f a c 1 _ . . 
In— = / dx——-+j 7 e^3{oit,Oi b ,a c ), 

(J, Ja 3 P3{X) 

in^= r dx 6 -^ + r dx 6 -^ + r + r 6 dx s -^ 

m i3 J a 3 fo(x) Ja c Pi(x) Ja b (3 b (x) J a t (3q{x) 

+ G^z{oi U a b ,a c ) . (41) 
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The dependence of our running couplings on mass logarithms can be deter- 
mined order by order from perturbative results for the Callan-Symanzik and 
matching functions. An example is the LO asymptotic behaviour of a t for 
simultaneous t, b, c decoupling: 

a t ~ 67r/(231n(m t //i) + 2\n(m b /fl) + 21n(m c //i)) . (42) 

Generally, this coupling runs with a mixture of / = 3, 4, and 5 active flavours. 
The sequential limit (4) is a special case, where the dominant term in the de- 
nominator exhibits the five-flavour coefficient 23 (33 — 2f for / = 5). Another 
extreme case is m t = m b with ln(m c / p) neglected, where the coefficient be- 
comes 25 (/ = 4). Similarly, we get the / = 3 coefficient 27 for m t = m b = m c . 

The multi-scale RG can be readily generalized to include logarithmic momen- 
tum scales between the mass logarithms in Eq. (2). The main new feature is 
that the invariant masses given by Eqs. (35) and (40) are used to construct 
scale-invariant logarithms. For example, let Q be a large momentum scale 
contained in the simultaneous limit 

\n(m b /p) < \n(Q/p) < \n(m t /p) -> large. (43) 

Then amplitudes are expanded as multinomials in a t , a b and «q, where cvq is 
defined by 

\n(Q/m b ) = f aQ dx/(3 5 (x) (44) 

Ja b 

or equivalently by 

rat 

Hm/Q) = / dx/fc(x) . (45) 

JOiQ 

This completes the formulation of the multi-scale RG. Its advantage over the 
standard method (1) is that to a given order, no logarithms are discarded, so 
uncontrolled approximations or outright contradictions (such as the example 
of anomaly cancellation below) are avoided. Each application of the standard 
method should be checked to see how well it approximates the multi-scale 
result. 

The need for a multi-scale approach is best illustrated by the classic exam- 
ple [7,13,14,15,16,17] of heavy- quark decoupling from the weak neutral axial- 
vector current 

J% = \{ E - E }(ffVM) 6 (46) 

Kq=u,c,t q=d,s,b) 

where the subscript 6 indicates an operator or amplitude in QCD with / = 6 
flavours. Diagrams for LO contributions in ghost-free gauges are shown in 
Fig. 1. 
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Fig. 1. Leading-order diagrams for heavy quarks h to decouple from light quarks 
I when the weak neutral current is coupled to a nucleon at low momentum 
transfer. The gluon propagators (indicated by wavy lines) are dressed with multiple 
one-loop self-energy insertions. 

Let us apply the standard procedure (1) to the six-flavour amplitude 

rg = (*vytf - h»isb) 6 . (47) 

The first step m t — > oo produces a current operator with scale-invariant 
renormalization [14,15]: 



lim r<g = rg = -((6 7M75 6) inv ) . (48) 



In the LO correction to this result, 

, (6) (5 ) _ 6a t I 



the Witten running coupling a* displays a /we-flavour dependence [7,14]: 

a t ^— — . (50) 

231n(m t //i) V ; 

This agrees with the rule of thumb that the number of active flavours is de- 
termined by the residual theory. 

The problem becomes more subtle when both t and b are to be decoupled 
such that anomaly cancellation occurs for m t = mj,. The sequential method 
(1) requires us to decouple the b quark from Eq. (49), 



r S LO = Seq 2^ ^ ~ ^ ( £(^ 759 ) ln v) 4 ^ 



but it is immediately evident that the method fails with respect to cancellation 
at m,b = m t because on, runs with four flavours: 

6n 

Qb ~ 251n(m^) " (52) 

Presumably this is why the sequential method was not pursued for this ex- 
ample in [7,14]; instead, a five-flavour dependence was proposed for both the 
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t and b contributions: 



r S 1 (s) 2 {r^m " u^m Xi^'U, ■ (53) 

Certainly this formula vanishes for m t = m fe , but the five-flavour dependence 
assumed for both scales is not possible: the residual theory has only four 
flavours. 

The multi-scale RG resolves these problems very simply. Eq. (51) is a double 
asymptotic series of the form (21) derived for the strictly sequential limit (1). 
One cannot apply it directly to the case m t = because important LO 
logarithms have been discarded. To restore all logarithms, replace Witten's 
running couplings a t , <S& by our simultaneous couplings a t , a& and so obtain a 
simultaneous expansion of the form (22): 

r 2»^(*-*)(£(^L) 4 - < 54 > 

Since 6 a t = a b for m t = m b , anomaly cancellation is automatic; at the same 
time, this result includes the sequential formula (51) as a special case. These 
properties are especially evident in the asymptotic result replacing (53): 

(tMl)„,) 4 . (55) 

x q=u ' 4 

A direct analysis of the diagrams in Fig. 1 gives the same answer [16]. 

The next order in the simultaneous expansion can be deduced from the NLO 
correction to the sequential result (51) found in [15]. Making the replacement 
ott,b a t,b, we find: 

U6) 6 / xf 125663 6167 1/w \ \ frn , 

The factor a b — a t ensures that anomaly cancellation occurs automatically at 
m t = m b . 

In practice, the multi-scale RG will be important in determining momentum 
dependence between heavy-quark thresholds, e.g. in deep inelastic scattering. 
It will allow the whole range of momenta, from one threshold to the next, to 
be described by a single expansion calculable (in principle) to any order. 

When a threshold is crossed, the multi-scale asymptotic expansion changes 
because the relevant operator-product and heavy-quark expansions differ. 
Consequently there are slope discontinuities at thresholds, as for other mass- 
independent methods [18]. The standard prescription for the removal of dis- 



3b 

sim 23 



r (6) 

1 ..c ~ 



251n^r 
u J 



231n^ + 21n^ 



u 



a 



10 



continuities is the method of effective charges [5,19,20,21], but this has been set 
up only for the single-scale RG. Formulating a natural multi-scale extension 
of this method presents an interesting challenge for the future. 

The result of this letter is a substantial generalization of the RG which solves 
problems involving logarithms for more than one large scale. Amplitudes are 
expanded as multinomials of running couplings, one for each large scale but 
depending (in general) on all of these scales; this replaces the familiar per- 
turbation series in one running coupling for a single-scale RG. All multi-scale 
results to a given order (degree in the multinomial expansion) can be readily 
deduced from single-scale calculations to the same order. Our method avoids 
inconsistencies by retaining all important logarithms to a given order. 
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